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Abstract 

A non-destructive, gamma-spectrometric method for verifying the plutonium con- 
tent of Pu-Be neutron sources has been developed. It is also shown that the isotopic 
composition and the age of plutonium (Pu) can be determined in the intensive neu- 
tron field of these sources by the "Multi-Group Analysis" method. Gamma spectra 
were taken in the far-field of the sample, which was assumed to be cylindrical. The 
isotopic composition and the age of Pu were determined using a commercial imple- 
mentation of the Multi-Group Analysis algorithm. The Pu content of the sources 
was evaluated from the count rates of the gamma-peaks of 239 Pu, relying on the 
assumption that the gamma-rays are coming to the detector parallel to each other. 
The determination of the specific neutron yields and the problem of neutron damage 
to the detector are also discussed. 

Key words: NDA, gamma-spectrum, Pu-Be source, Pu content, Pu isotopic 
composition, Age dating 
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1 Introduction 



In recent years Pu-Be sources have been intensively studied [1,2,3] because 
they represent a nuclear safety and safeguards issue. A large number of such 
sources (around 200) - mostly out of use - are stored in Hungary (and this is 
the case in several neighbouring countries as well). The Pu-content of these 
sources is to be accounted for, regularly reported to and inspected by the 
International Atomic Energy Agency, and from 2004 on, they are also subject 
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to EURATOM safeguards. Because the late Soviet supplier did not declare the 
Pu-content of the sources, the values for the Pu content presently on record are 
based on rough overestimates assuming that the sources contain pure 239 Pu. 

Using the Multi-Group Analysis algorithm [4,5] to analyze the gamma spectra 
of some Pu-Be sources, however, has shown that the 239 Pu component varies 
in the range of 75-95%. It is known that the specific neutron yields can be 
very different for different Pu isotopes, for example, ~ 10 4 n/sg for 239 Pu, and 
~ 10 7 n/sg for 238 Pu as calculated [2] from specific alpha yield [6,7]. Conse- 
quently, the total Pu content should be corrected by using the following data: 
1. the measured Pu isotopic composition, 2. neutron output from declared or 
measured data and 3. specific yield values. The calculating procedure which 
uses these corrections we call "combined method" or "CBM" and its results 
were reported in Refs. [1,2,3]. In case of lower 239 Pu isotope ratio (about 70- 
80%) the total Pu contents are substantially lower than the values still kept in 
the files for safeguards. A typical example is a source with 24 g of Pu found by 
the CBM instead of the value of 178 g on file. However, since specific neutron 
yields are not isotopic physical parameters, they are strongly dependent on 
the technology applied by the manufacturer. For example, the value for the 
specific neutron yield of 239 Pu may vary in the range (1.2 — 6.2) x 10 4 n/sg, 
[7]. Therefore, the combined method may contain systematical errors because 
of the lack of the precise knowledge of the Pu specific neutron yields. The Pu 
contents estimated by the CBM should be considered as upper limits for the 
Pu content [2]. Similarly, the lower limits can be estimated if the maximum 
specific neutron yields are used. Although these estimates are subject to the 
high uncertainty associated with the specific neutron yields needed for the cal- 
culation of the Pu-content, the CBM gave an indication that the Pu contents 
of Pu-Be neutron sources must be carefully investigated. 

In parallel with the gamma-spectrometric measurement of the isotopic compo- 
sition, a high-resolution gamma-spectrometric method which does not require 
the knowledge of the specific neutron yields was developed for the experimen- 
tal determination of the plutonium content of the Pu-Be sources. Further- 
more, combining the results of the gamma-spectrometric measurement with 
the measurement of the total neutron output the specific neutron yields of the 
sources can also be determined, thus providing additional information about 
the sources. 

During developing this topic, another approach using neutron coincidence 
counting technique was considered for the estimation of the Pu content. A 
neutron coincidence counter was built and a correlation was obtained between 
the Pu content and the ratio of the number of real coincidences (R) to the to- 
tal number of counts (T). This method will be called here the "R/T" method. 
An important advantage of this method is that there is no need for using 
the Pu isotopic composition measured by gamma spectrometry. However, to 
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develop the R/T method standard Pu-Be neutron sources of known isotopic 
composition and Pu content are needed for getting an accurate calibration 
curve. Unfortunately, there are no such sources in our laboratory and, as far 
as we know, there are none in the laboratories with which we made contact. 
The gamma-spectrometric method described here made it possible to prepare 
a set of calibrating sources which can be used for developing the R/T method. 
The results of this research can help to develop the R/T method and to mon- 
itor the specific neutron yield of the sources. 

This paper presents a procedure for evaluating the Pu contents of Pu-Be 
neutron sources based purely on high-resolution gamma-spectrometry. Pu- 
Be sources with neutron output from 10 4 n/s to 10 7 n/s were assayed. The 
samples were set up in the far field of the detector, so that the gamma-rays 
may be assumed to be coming to the detector parallel to its axis. As the Pu- 
Be sources are usually encapsulated in steel cylinders of known wall thickness, 
the correction factor for self-absorption can be appropriately calculated by a 
numerical model for the far-field assay of a cylindrical sample. 

Note that the Pu content can also be derived from the count rates of the 
gamma-rays of 239 Pu by a method called the "infinite energy method" [8]. For 
Pu-Be sources, however, the uncertainty of that method may reach up to 20-25 
%. If, however, the self-absorption correction can be accuretely calcualted, as 
is the case with the samples assayed in this work, the precision of determining 
the Pu content by the procedure described in this paper is expected to be 
better than the precision of the infinite energy method. 



2 Determination of the isotopic composition and of the age of plu- 
tonium 

Spectra of Pu-Be neutron sources were taken in the 0-600 keV energy region 
by a 2000 mm 2 planar high-purity germanium detector (Canberra GL2020R) 
at 20-200 cm source-to-detector distance (see Fig. 1). The parameters of the 
sources and of the measurements are given in Table 1. The multi-channel 
analyzer was set up at 8196 channels to get the value of 0.075 keV per channel 
as a requirement of applying the Multi-Group Analysis computer code (MGA) 
[4]. The plutonium isotopic composition was determined by the commercial 
implementation of MGA, MGA++ [5]. The age is estimated by a software on 
the basis of the daughter /parent activity ratio of 241 Am/ 241 Pu. The isotopic 
composition and the age of 8 Pu-Be sources, which are mainly the sources used 
in a previous paper [2], are shown in Table 1. The values obtained by MGA 
for the isotopic composition are in agreement with those obtained in Ref. 
[2] by analyzing gamma- and x-rays below 300 keV. The age calculated by 
MGA agrees with the one which can be obtained from the date of production 
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available in the records. 



HPGe detector 



a) 



Pu-Be source 
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Container 
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Fig. 1. The measurement setup 

The isotopic fractions of 238 Pu, 240 Pu, 241 Pu, 242 Pu and 241 Am are plotted 
versus 239 Pu isotopie ratio on Fig. 2. It can be seen that there is a correlation 
between them. This correlation can be explained by the production of Pu 
isotopes in a nuclear power reactor. 



3 Determination of the Pu content 



3. 1 Calculating procedure 



A typical gamma spectrum of a Pu-Be source is shown in Fig. 3, where the 
129, 203, 345 and 413 keV lines of 239 Pu can be clearly seen. 

The correlation of the count rate of a gamma-ray with energy E, C E , and 
total Pu content, M Pu) can be expressed as 

C E (M Pu ) = M Pu f 239 G E E ^r , (1) 
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Fig. 2. Correlation between the 239 Pu isotope ratio and isotope ratio of other Pu 
isotopes and 241 Am 

where ^239 is the 239 Pu isotopic ratio, Ge is the specific gamma yield (gamma/gs) 
[7], Oe is the efficiency of detector measured by standard sources and Fe is 
a correction factor for absorption. F E is the product of a factor accounting 
for self- absorption, Fie, and a factor accounting for absorption in the steel 
wall of the container, F 2 e- Assuming that the gamma-rays come parallel to 
the detector axis, Fie can be calculated using the far-field approximation for 
a cylindrical sample viewed along a diameter [7,9]: 



1 



'IE 



2 h(v 1E D) - L^eD) 



(2) 



while 



F 2E = exp(fi 1CE d c ) , 



(3) 



where 



Hie is the linear attenuation coefficient of the Pu-Be material, 

D is the diameter of the sample, 

Ii is modified Bessel function of order 1 [10], 
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Fig. 3. Gamma spectrum of a Pu-Be source taken by a planar high-purity germanium 
detector 

• L\ is modified Struve function of order 1 [10], 

• Vice is the linear attenuation coefficient of the steel container taken from 
Ref. [11] and 

• dc is the thickness of the container wall. 

Furthermore, can be calculated from the mass attenuation coefficient 
(cm 2 /g) of Pu, fi PuE , and Be, p BeE , [11] as: 



Vie = PVe 



239 



239 + 9n 



VPuE + 



9n 



239 + 9n 



VBeE 



(4) 



where p is the density of the sample, n is the ratio of the number of Be atoms to 
the number of Pu atoms (PuBe„), and the numbers 239 and 9 are the atomic 
mass numbers of Pu and Be, respectively. In general, the Pu-Be material 
probably does not completely fill the inner volume of container. Nevertheless, 
it is expected to take on the cylindrical shape of the container, so the diameter 
D of the Pu-Be material can be expressed as: 



D(p,R,n,M Pu ) 



AM Pu 239 + 9n\ 1 /3 

x 

npR 239 / 



(5) 
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where R = H/D is the ratio of height and diameter of the Pu-Be material 
(see Fig. 1). 

The calculated values of F E = F 1E x F 2E , 1/F E and C E for the 129, 203, 345, 
375 and 413 keV lines are plotted versus plutonium mass on figures 4, 5 and 6, 
respectively, in the range below 120 g. The values of the parameters were set 
as n = 13, p = 3.7 g/cm 3 for a "usual" Pu-Be material, and it was assumed 
that R — 1, dc — 4 mm, and ^239 = 1 for a "model" source. Furthermore, 
L = 50 cm was taken as a typical measurement distance, while l\ and L\ 
were taken from Ref. [10]. Using these values for the parameters, the following 
observations can be derived from formula (1) and Figs. 4-6: 




20 40 60 80 100 120 140 

239 Pu mass (g) 

Fig. 4. Dependence of the correction factor for absorption, F E , on the mass of 239 Pu, 
assuming a cylindrical sample. 

The value of F is relatively small for the 345, 375 and 413 keV lines and it 
varies slowly between 1.5 and 3 for plutonium masses in the range 0-120 g. 
Consequently, about 30-70 % of the gama photons at these energies leaves the 
source. On the other hand, for the 129 keV and 203 keV lines, the value of F 
is relatively large. Furthermore, it rises very quickly (from 3 to 80 for 129 keV 
and from 2 to 10 for 203 keV for plutonium masses in the range 0-120 g) and 
only a small percentage of these gamma rays escapes from the source (1-30 % 
for 129 keV and 10-50 % for 203 keV). Note finally, that the count rate of the 
129, 375 and 413 keV lines is much higher than that of the 203 and 345 keV 
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Fig. 5. Calculated values of I/Fe versus total plutonium mass, Mp u , for five 239 Pu 
gamma-rays. 

lines. The first group will be called "intensive group". 

For determining the plutonium content the use of gamma peaks of relatively 
high intensity is preferred, the count rate of which can be measured with a 
satisfactory statistical error within a reasonable measurement time. Namely, 
since neutron damage to the detector crystal is to be reduced to the minimum, 
as short a measurement time as possible should be chosen. Furthermore, the 
count rate of the chosen gamma peaks should not be very sensitive to the 
geometrical dimensions of the assayed sample, because the latter are only 
known with a relatively large uncertainty. The gamma peaks at 375 and 413 
keV satisfy both these requirements. The count rate of the peak at 129 keV 
is also high, but it strongly depends on the geometrical dimensions of the 
sample, so in the present work it has not been used for the quantitative assay 
of the samples. Note, however, that exactly this property of the 129 keV line 
can be used for verifying the parameters of the source, i.e., for estimating the 
ratio R = H/D. This aspect will be reported in a subsequent paper. 

Equation (1) can be rewritten in the form 

M Pu = - %- F E = f(C E , / 239 , G E , O e , P) , (6) 



8 




10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 

Total plutonium mass M pu (g) 

Fig. 6. Calculated count rate of the peaks of 239 Pu at 50 cm distance from the 
source. 

where P is the set {n, p, R, d c } of parameters of the source. Using formula (6) 
M Pu will be derived from the count rates Ce by taking an adopted value for 
the parameters in the set P. It means that for each source there will be more 
values for Mp u , calculated from the count rate of gamma peaks at different 
energies. In the present case we will only consider the values of Mp u which 
can be obtained from the count rates of the gamma peaks at 375 and 413 keV. 
The two values should agree with each other and the final value of Mp u will 
be calculated as the weighted average of these two, for each source separately. 



3.2 Analyzing the sources of errors 

Although, in general, each Pu-Be neutron source has a certification document 
providing the parameters of the container (see Table 1), the values in the 
parameter set P are not always exactly known. The error of Mp u , AM, caused 
by the uncertainties of these parameters can be estimated as 



W . (7) 



AM 



AC f 



F E = 



J239GeOe Jiz'aGeOe ^ 
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where Pk denotes the fc-th element in the parameter set P. The uncertainties 
associated with systematical errors of Ge, Oe and Fe are ignored in the above 
equation, because in this section we wish to concentrate only on the sources 
of error connected to the inappropriate knowledge of the parameter set P. 

The values of Ce calculated according to equation (1) are plotted on Figs. 7 
and 8 for the 129 and 413 keV lines as functions of the parameters n, p, R 
and dc, for two typical Mp u values of 0.02 g and 60 g. The values of Ce for 
413 kev vary very slowly with n and p, and a bit more quickly with R and 
dc- This means that the values of dCE/dPk m equation (7) are small in the 
case of 413 keV. The same is true for 375 keV as well. Therefore, the total Pu 
content can be estimated from the count rate of the 375 and 413 keV peaks 
with high accuracy even if one uses the uncertain values of the parameters 
given in the certificate of the Pu-Be source. A rough estimation based on 
the calculation model presented above and illustrated in Figs. 7 and 8, using 
typical parameters of the sources, yields that the relative error of Mp u caused 
by the uncertainties of these parameters is below 5 %. 

The value of C E for the 129 keV peak is more sensitive to the changes of R 
and dc, thus the uncertainties of R and dc cause a larger error in the value of 
Mp u calculated using the count rate of the 129 keV peak. Therefore, as already 
remarked, in this work the 129 keV peak will not be used for estimating the 
plutonium content. 

3.3 Determination of the Pu content 

The parameters of the sources, the measurement set up, the count rates and 
the Pu contents calculated using formula (6) are presented in Tables 1-5. The 
final Pu contents (Table 5, column 6) are the averages of the Pu contents 
estimated from the count rates of the 375 and 413 keV lines. 



4 Results and discussion 

The Pu contents obtained by the described method are smaller by about 30- 
50% than the values obtaied previously for the same Pu-Be sources by the 
method which uses the isotopic composition and the total neutron output 
(the "combined method" or CBM [2]) for calculating the Pu content. This 
implies that the values of the specific neutron yield are larger than the ones 
used in the CBM, by about 30-50%. If, in addition, the number of neutrons 
from (ct,n) reactions is estimated, for example, by measuring the 4438 keV 
gamma line [12], then the specific neutron yield can also be determined, using 
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Fig. 7. Calculated count rate of the 129 and 413 keV peaks of 239 Pu at 50 cm 
distance from the source as a function of various parameters of the source with 
total plutonium mass Mp u = 0.02 g. 
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Fig. 8. Calculated count rate of the 129 and 413 keV peaks of 239 Pu at 50 cm 
distance from the source as a function of various parameters of the source with 
total plutonium mass Mp u = 60 g. 
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the results of the present work. 



The calibration curve for the M — R/T correlation of the "R/T" method 
mentioned in the introduction of this paper becomes more accurate if one 
uses the total Pu contents obtained by the present method instead of the ones 
obtained by the CBM. 

In each measurement, the detector was exposed to an integral neutron flux of 
about 10 4 — 10 5 n/cm 2 to get a statistical error below 5% for the peaks in the 
"intensive group", IG. The measurement time should be much longer if the 
statistical error of all five relevant gamma-rays of 239 Pu is to be bellow 5%. 
In the described measuremts the degree of neutron damage to the detector is 
expected to be negligible, because in the present measurements the integral 
neutron flux was much smaller than 10 9 n/cm 2 , which is the integral flux that 
causes a noticeable damage to the detector, as estimated by some authors 
[13,14]. 

The plutonium content of the sources can also be estimated by the the so-called 
"infinite energy method" [8] in the following way. The value of m(l/E) = 
C e / {fvsvG eO e) is plotted as a function of 1/E for five gamma-rays (129, 203, 
345, 375 and 413 keV). The curve fitted to these points is then extrapolated 
to 1/E = 0, and the intercept, m(l/E = 0) yields the plutonium content 
of the source Mp u . Comparing m(l/E) to formula (1), it can be seen that 
m(l/E) = Mp u /Fe- Therefore the relation between the method presented in 
this paper and the infinite energy method is that in this work Fe was used 
to correct for self-absorption, instead of using the extrapolation of the infinite 
energy method. Calculating F E requires the knowledge of certain parameters of 
the sources, such as, e.g., the dimensions of the source, the type of the container 
etc., whereas the infinite energy method does not. However, if the required 
data are available, the method presented here gives a more accurate result for 
the plutonium content than the infinite energy method. In the present case 
these data were available from the certificates of the assayed sources, so their 
plutonium content was evaluated using the method described in this work. 
An exception is the source having the identifier 4.7-555 (Table 5) which was 
encased in a double container of unknown dimensions, so its plutonium content 
was estimated using the infinite energy method. 

Alternatively, if the source dimensions are not exactly known, the 129 keV line 
can be used for monitoring the parameters of the sources. Presently a method 
based on using the 129 keV line is being developed in order to determine the 
plutonium content of completely unknown Pu-Be sources with high precision. 
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Note added in the revised version 

After the original version of this paper has been submitted, the results of the 
gamma-spectrometric plutonium-mass measurements of several Pu-Be sources 
were confirmed by calorimetric measurements [15]. The comparison of the two 
methods will be dealt with in a subsequent paper. 
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Tables 



Table 1 

Parameters of the assayed Pu-Be sources taken from the certification documents and the parameters of measurements. The first numbc 
in the source identifier is its nominal Pu content in grams. (*Put in an aluminum container.) 





JNeutron 


Outer dimensions 


Inner 


dimensions 








Identifyer 


output 


of the 


container 


of the 


container 




Measurement parameters 




(n/s) 


(cm) 




(cm) 














Diameter 


Height 


Diameter, D Height, H 


I 

(cm) 


Meas. time 
(s) 


JNeutron flux 
(xlO 4 n/cm 2 ) 


0.2-480 


1.1 x 10 4 


1.0 ±0.02 


1.9 ±0.05 






20 


680.0 


0.15 


2-479 


1.1 x 10 5 


1.5 ±0.02 


2.4 ±0.02 






50 


1352.85 


0.47 


1.9-611a 


1.16 x 10 5 


2.0 


3.0 


1.3 


1.95 


50 


1198.9 


0.44 


4.7-555* 


2.9 x 10 5 


2.0 


3.0 


1.3 


1.95 


50 


695.5 


0.59 


4-407 


2.68 x 10 5 


2.0 


3.0 


1.3 


1.95 


50 


974.0 


1.02 


37-425 


2.26 x 10 6 


2.4 ±0.02 


3.0 


1.6 


1.60 


200 


3326 


5.71 


178-461 


1.1 x 10 7 


3.5 


4.5 


2.7 


2.70 


200 


606.5 


1.33 


85-701 


5.27 x 10 6 


3.5 


4.2 


2.9 


3.15 


200 


648.8 


0.68 



Table 2 



Isotopic composition and age of the sources determined by the MGA code from the gamma spectra taken in the region 0-600 keV. 



Identifyer 






Isotope ratio (%) 






Age (years) 


Date on 
certificate 


Measure- 
ment date 




238p u 


239p u 


240 Pu 


241 Pu 


242p u 


241 Am 








0.2-480 


1±3.5% 


83±1% 


14±5% 


0.71±4.4% 


1.19±10% 


2.2±3.5% 


28.4±0.6 


1982.10 


2002.07.25 


2.479 


1±3% 


79.9±1% 


17±4% 


0.73±3.5% 


1.38±10% 


2.27±3% 


29.7±0.5 


1997.10 


2002.07.25 


1.9-611a 


0.004±50% 


95.6±0.2% 


4.3±4% 


0.046±5.5% 


0.014±10% 


0.23±2% 


37.6±1 


1965.06.08 


2002.06.25 


4.7-555 


0.0035±75% 


96±0.2% 


4.15±4.5% 


0.05±6.5% 


0.012±10% 


0.22±2.5% 


35.5±1.3 


1967.03.16 


2002.02.23 


4-407 


0.001±200% 


95±0.2% 


4.73±3.5% 


0.04±6.5% 


0.015±10% 


0.23±2% 


39±1.3 


1965.01.30 


2002.11.19 


37-425 


1.06±3% 


77±1% 


19.7±4% 


0.7±4% 


1.6±10% 


2.8±3% 


33.5±0.5 


1973.10.29 


2004.12.09 


178-461 


1.3±4% 


76±1.5% 


20±4% 


0.8±6% 


1.9±10% 


2.9±4% 


31.5±0.5 


1976.08.26 


2002.04.28 


85-701 


0.0076±56% 


95±0.3% 


5±6% 


0.044±12% 


0.017±10% 


0.22±3.5% 


37±1.5 


1966.08.20 


2002.07.25 



Table 3 

Parameters used for calculations 



Gamma-energy 


129 keV 


203 keV 


345 keV 


375 keV 


413 keV 


HE (cm 2 /g) [11] 


2.05 


0.9 


0.3 


0.27 


0.205 


W (Fe) (cm- 1 ) [11] 


1.85 


1.15 


0.81 


0.77 


0.73 


G E (10 4 photons/gs) [7] 


14.361 


1.2847 


1.2829 


3.6014 


3.6 


Oe x 10 4 (measured at 50 cm) 


4.6 


2.8 


1.36 


1.21 


1.07 



Table 4 

Count rate of 239 Pu gamma-rays 



Identifier Count rate (cps) 





129 keV 


203 keV 


345 keV 


375 keV 


413 keV 


0.2-480 


2.25±0.12 


0.191±0.059 


0.159±0.041 


0.31±0.059 


0.307±0.059 


2-479 


1.83±0.11 


0.239±0.037 


0.174±0.032 


0.448±0.04 


0.427±0.044 


1.9-611a 


5.88±0.13 


1.053±0.083 


0.799±0.05 


2.08±0.06 


1.92±0.063 


4.7-555 


13.63±0.36 


2.21±0.21 


2.21±0.15 


5.19±0.18 


4.66±0.15 


4-407 


14.55±0.21 


2.36±0.12 


2.15±0.08 


5.75±0.13 


5.21±0.12 


37-425 


0.85±0.05 


0.15±0.03 


0.14±0.015 


0.38±0.02 


0.375±0.017 


178-461 


2.67±0.26 


0.53±0.16 


0.63±0.13 


1.53±0.15 


1.52±0.10 


85-701 


5.79±0.29 


1.15±0.15 


1.35. ±0.12 


3.96±0.15 


4.02±0.15 
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Table 5 

Total plutonium content estimated from the count rate of the 375 and 413 keV lines. 
Identifier Parameters Total plutonium content (g) 





R 


dc (cm) 


375 keV 


413 keV 


Average 


0.2-480 


1.9 


0.35 


0.019±0.004 


0.021±0.004 


0.020±0.004 


2-479 


1.6 


0.35 


0.20±0.02 


0.20±0.02 


0.205±0.04 


1.9-611a 


1.5 


0.35 


0.85±0.03 


0.83±0.03 


0.84±0.03 


4.7-555 


unknown 


infinite ener 


gy method 


3±0.4 


4-407 


1.5 


0.35 


2.65±0.07 


2.47±0.06 


2.55±.05 


37-425 


1 


0.4 


4±0.25 


4±0.2 


4.0±0.2 


178-461 


1 


0.4 


19.3±2.3 


19.5±2.2 


19.5±1.6 


85-701 


1.08 


0.3 


49.5±2.5 


49.1±2.3 


49.2±1.5 
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